Spin Dynamics of a J\ — J2 antiferromagnet and its implications for iron pnictides 
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Motivated by the recent observation of antiferromagnetic correlations in the paramagnetic phase 
of iron pnictides, we study the finite-temperature spin dynamics of a two-dimensional Ji — J-z 
antiferromagnet. We consider the paramagnetic phase in the regime of a (7r,0) collinear ground 
state, using the modified spin wave theory. Below the mean-field Ising transition temperature, 
we identify short-range anisotropic antiferromagnetic correlations. We show that the dynamical 
structure factor <S(q, to) contains elliptic features in the momentum space, and determine its variation 
with temperature and energy. Implications for the spin-dynamical experiments in the iron pnictides 
are discussed. 



I. INTRODUCTION 

High-temperature superconductivity in the iron 
pnictides^ arises by doping antiferromagnetic parent 
compounds^. Hence, the strength of the electronic cor- 
relations, the nature of magnetism, and the relationship 
between magnetic excitations and the superconductivity 
are important issues for understanding the emergence of 
high temperature superconductivity in these materials. 
In the parent iron pnictides, the Neel transition into a 
(-7T, 0) antiferromagnet is either preceded by or concomi- 
tant with a tetragonal-to-orthorhombic structural tran- 
sition. The (-7T, 0) magnetic order by itself can be under- 
stood either by invoking a local moment J1 — J2 modell- 
er an itinerant model with nearly nested electron and 
hole pocket*^— . 

The experimentally observed "bad metal" behavior, 
the Drude- weight suppressio n 15 ' 16 and the temperature- 
induced spectral- weight transfer— ~— place these materi- 
als near to a Mott transitio n 4 ! 9 ' 19 ! 20 ; a Mott insulator can 
emerge when the iron square lattice either expands^ or 
contains ordered vacancies 2 ^. In a metallic system close 
to a Mott transition, quasi-local moments are expected to 
arise; this picture is further supported by the experimen- 
tal observation of zone boundary spin wave excitations 
in the magnetically ordered state at low temperatures— 
The inelastic neutron scattering experiments demon- 
strated the need for an anisotropic J\ — J2 model with 
Jix 7^ Jiy, which may reflect an orbital ordering 2 ^ - — 
while pointing to the relevance of magnetic frustration 
from the extracted ratio {J\ x + Ji y )/2J 2 ~ 1—. There- 
fore, results in the tetragonal, paramagnetic phase of the 
parent compounds are of great importance for under- 
standing the relevance of an isotropic J\ — Ji model as 
well as the strength of the underlying magnetic frustra- 
tion. Recent inelastic neutron scattering measurements 
of Diallo et al£l on the tetragonal, paramagnetic phase 
of CaFe2As2 represent a first step in this direction. Even 
above the concomitant first-order structural and Neel 
transition temperature, they have observed anisotropic 
spin dynamics around the (jr, 0) wave vector, and the in- 
ferred ratio J1/J2 ~ 0.55 is similar to that of the ordered 
phase. 



Motivated by these experimental results we study the 
spin dynamics of a two-dimensional J\ — J2 antiferromag- 
net. While theoretical studies exist on the order- from- 
disorder phenomenon and phase diagram of the J\ — J2 
mode l 28 ! 29 , the spin dynamics in the paramagnetic phase 
of the model in the (tt, 0) collinear regime has not yet 
been systematically studied. We carry out the calcula- 
tions using a modified spin wave theory^, which incorpo- 
rates the 1/S corrections that are important for captur- 
ing the order-from-disorder phenomenon and the associ- 
ated dynamical properties. We discuss the implications 
of our results for the iron pnictides, including the role of 
itinerant electrons. 

Our paper is organized as follows. In Sec. II we intro- 
duce the relevant J\ — J2 model and describe the modified 
spin wave theory calculations. In Sec. Ill we analyze the 
excitation spectrum obtained from modified spin wave 
theory, and associated behavior of the spin-spin correla- 
tion length. In Sec. IV we analyze the dynamic structure 
factor calculated by using the modified spin wave theory 
results. In Sec. V wc consider the fluctuation effects due 
to itinerant electrons within a Ginzburg-Landau frame- 
work. In Sec. VI we describe the relation between our 
theoretical results and the experimental data obtained 
in the paramagnetic phase of iron pnictides. We provide 
a summary of our work in Sec. VII. The technical de- 
tails of fitting the experimental data and consideration of 
inter-planar exchange coupling using modified spin wave 
theory are respectively relegated to Appendix A and Ap- 
pendix B. 

II. MODEL AND MODIFIED SPIN WAVE 
THEORY 

The model is defined by the Hamiltonian 

h = j 1 ]Ts i -s i + J 2 Si ' s J> M 

m <«» 

where J\ and J2 respectively denote the antiferromag- 
netic exchange couplings between spins located in the 
nearest ((ij)) and next-nearest neighbor (((ij)}) sites on 
a square lattice. Classically, for J2/J1 > a c = 0.5, the 
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FIG. 1: (Color online) The temperature dependence of the 
mean-field parameters, for 5 = 1 and J1/J2 = 0.8. The 
decoupled Neel sublattices are illustrated in the upper right 
corner, which also defines the angle <j). 



izes their dispersion; in this respect, it is similar to the 
large- TV Schwinger boson mean-field theory 3 ^. Following 
Takahash i 30 ' 32 , we express the Hamiltonian, Eq. (1), in 
terms of the DM bosons in momentum space. The proce- 
dure is to minimize the free energy T = (H) — TS under 
the constraint of zero magnetization, (S—alai) = 0, with 
respect to variational parameters which enter J- . These 
are the boson dispersion ek, the angle <f> and the Bogoli- 
ubov angle 6^. The latter enters in a Bogoliubov trans- 
formation that mixes the operators of the two interpene- 
trating Neel sublattices and renders the nonzero temper- 
ature density matrix diagonal 3 ^. The equal time correla- 
tors (Si -Sj) can be written in terms of expectation values 
like {a\a,j). Therefore, we define ferromagnetic and an- 
tiferromagnetic bond correlations fij — (a\a,j) — (<Zi<zJ) 

and gij = (aiaj) = (ajaj). The explicit expressions for 
the bond correlations arc given by 



fij = cosh26» k (n k + -) exp(-ik • r y ) (2) 



lattice decouples into two independently Neel ordered, 
interpenetrating lattices, and the angle between the 
staggered magnetizations of these two sublattices, as il- 
lustrated in the Fig. [1] inset is arbitrary. An ordcr-from- 
disorder transition at temperature T a breaks the four- 
fold rotational symmetry of the square lattice down to 
a twofold rotational symmetry of the rectangular lat- 
tice, and (f) = 0, 7t emerge as degenerate ground states 
at T = 02£. Since quantum fluctuations make a c > 0.5, 
for definiteness we will focus on J2/J1 > 1. 

We define a local spin quantization axis along the clas- 
sical ordering direction at each site (£lf), as illustrated 
in an inset to Fig. [TJ We then introduce the corre- 
sponding Dyson-Maleev (DM) boson representation for 
the spin operators at each site: Si ■ ilf = S — a, t 
well as S+ = V2S(l - a\ai/2S)ai and S~ = \f2Sa % 
The modified spin wave theory 3 ^ treats the self-energy 
of the a-bosons as a static quantity, which rcnormal- 



as 

t 



5 « = E siim2 #k("k + ^) exp(-ik • ry), (3) 



where n k = [cxp(e k /T) — 1] 1 is the Bose occupation 
factor. 

In terms of fij and gij the equal time spin correlator 
(Si ■ Sj) can be expressed as 



(Si • S 3 -) = cos 2 ?f 

„• 2 4>ij 



S + \- /(0) + fij 



S+^-f(0) + ga 



(4) 



where <f>ij = (j), ir — (f>, n, for horizontal, vertical and 
diagonal bonds respectively (see Fig. [TJ. Using the ex- 
pression for (Si • Sj) for different bonds, the total energy 
can be written as 



Si— ±£ 



cos 2 £ ( S + £ - /(0) + U ) - sin 2 Zl S +l-f(0) + g x 



S+^-f(0)+9y 



JiN 



E 

S 2 =±y 



Sin — X 
2 



J 2 N 



E 



S3=±x±y 



S+^-f(0)+g x+y 



(5) 



Notice that the expression for total energy only contains 
the nearest and next-nearest neighbor bond correlation 
parameters f x , f y , g x , g y and g x +y The constraint of 
zero magnetization, appropriate for T > T/v (for the two 
dimensional problem T/v = 0), is enforced by the La- 



grange multiplier fi. Minimizing E — TS — /^/(0) with 
respect to e k , <f>, # k , we obtain tanh20 k = A k /_B k , 
^ = V B l ~ K and sin0 (/ 2 + g\ - / 2 - .g 2 ) = 0, 
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where 



A k = 2Ji ( sin 2 -g x C x . k + cos 2 -g y C y .y 



-4J2 g x +yC x +yji 



sin 2 ^(9x-fy) + cos 2 ^(g y - 



+2Ji |^cos 2 -f x C X M + sin 2 -fyCy^J + 4J 2 .g^+y - Mi 

(7) 

and we have introduced the form factors C x ^ = cosk x a, 
C y .k = cos k y a, and C x + y y = cosk x acosk y a. Now using 
tanh2#k = A^/Bu in Eq. 02 we obtain the following set 
of self-consistent equations 



fa = 4f$3~ ("k + ^)c Q:k , a = x,j/ 



k 



(8) 



7a = 77 ^2 (" k + ^ ) 



, rikH — 
iV V £ k V 2 

k 



(10) 



We identify two important temperature scales To and 
T a Q such that To > Tjo, by solving the self-consistent 
equations. The temperature T = J 2 (S + l/2)[log(l/5 + 
marks the onset of the largest bond correlation 
g x + y , while T CT o marks the onset of nearest- neighbor bond 
correlations. For T > Tq, all the bond correlations vanish 
and we have decoupled local moment behavior. The first- 
order transition from the correlated to decoupled moment 
state at To is an artifact of the mean-field theory^. In the 
temperature range T a Q < T < To, the sublatticc angle 
(j) remains arbitrary, and the system has C^ v rotational 
symmetry. For T < Tro there are two degenerate solu- 
tions = 7r, with g y = fx = 0,g x ^ 0, f y ^ 0,g x ^ f y , 
and = 0, with x f-> y switching. An Ising order param- 
eter, which is defined classically as a = il\ ■ Q 2 = cos0, 
is modified to o oc 2 (cos 2 |(/| + g 2 y ) - sin 2 |(/ 2 + <? 2 )) , 

and becomes nonzero below T CT o- We identify this tem- 
perature as the mean-field "Ising transition" tempera- 
ture; fluctuations will reduce the actual transition to 
Tj < T CT o- In the following, we will focus on the state 
with </> = 7r. The spectrum is gapped at any nonzero 
temperature, but becomes gapless at T = giving rise to 
(it, 0) antifcrromagnetic order via a Bose condensation. 



III. LOW ENERGY SPECTRUM AND 
CORRELATION LENGTH 

The boson dispersion ek is shown in Fig.[5J For <fi = it, 
and T <C Tjo, the low energy physics is governed by the 
excitations in the vicinity of the ordering vector (it, 0) , 
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FIG. 2: (Color online)The dispersion £k along high symme- 
try directions in the paramagnetic Brillouin zone for differ- 
ent temperatures and 5 = 1, J1/J2 = 0.8. The curves 
from top to bottom viewed at the left end are for T/J2 = 
0.5,1.0,2.0,2.1,2.2. The plotted directions in the Brillouin 
zone are displayed in the upper right corner 



where the absolute minimum of the dispersion is located. 
Near (7r,0), the dispersion can be approximated by 



e k = [v 2 lx ^-k x ) 2 + v 2 ly kl + Aiy 

1 

Ai = [-p(8J 2 g x+y + 4Ji# x - /j)] 2 
v\ x = a(4J 2 g x+y + 2Jig x ) 

1 (AJ 2 g x+y + 2Jig x )(4J 2 g x + y - 2J\f y ) 



Vly 



+ 2Jlfyll 



(11) 

(12) 
(13) 



(14) 



Similarly in the vicinity of (0,7r), the excitation can be 
approximated as 



£ k = [vl x k 2 x +vl y (7T-ky) 2 +Al^ (15) 

A 2 = [(8J29x+y - 4Jif y - »)(4J l9x - 4/„ -//)]' (16) 
v 2x = a(4J 2 g x+y - 2Jxg x ) (17) 

v 2y = a 4J 2 g x+y (4:J 2 g x+y ~ 2Jig x ) + 2Jif y (4:J 2 g x+y 

1 

+2j l9x - aj^ - n)a&) 

At low temperatures T <C T CT o, the Lagrange multi- 
plier \x is exponentially small, and Ai <C A 2 . Therefore 
the spin-spin correlation length at low temperatures will 
be dominated by the smallest gap Ai = Texp[— Aj/T], 
where A ,7 = 2irp is the Josephson energy, with p = moV\ y 
being the stiffness and mo the staggered magnetization 
at T = 0. The velocity anisotropy yields two correlation 
lengths, £ x = Vi x /Ai and £ y = viy/At. 
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The low energy excitations around (ir, 0) can also be 
described in terms of an anisotropic 0(3) nonlinear sigma 
model. Ignoring the 1/5 corrections and weak temper- 
ature dependence of the bond parameters, we can take 
9x = fy = 9x+y = S, and obtain bare parameters of the 
sigma model xJjo = 4 (2J2 + Ji)a 2 , p x0 = (2J 2 + Ji)S 2 , 
and p y o = (2J2 — Ji)S 2 . The spatial anisotropy is 
captured by two direction dependent spin stiffness con- 
stants p x o and p y o, and x±o is the bare uniform trans- 
verse susceptibility. The spin wave velocities before 1/S 
corrections are given by V\ x = \J p x o/x±o> an d v iy = 
y/p y o/x±o- The temperature dependence of the gap is 
determined by the bare Josephson energy scale 



2ir Po = 4tt J 2 S 2 



/2 



(19) 



where po = ^PxoPyQ is the bare, geometric mean stiffness 
constant. For parameter values S = 1, and J1/J2 = 
0.8, we find Ajo = II.5J2. After solving the mean field 
equations, we obtain the Josephson energy scale 



Aj = = nm 



(4J29x+y + 2Jig x )(4:J2gx+y 



-2 Ji/„) + 2 Jrfyn 



(20) 



where m is the staggered magnetization at zero tempera- 
ture and captures the 1/S corrections to Aj. For S = 1, 
J1/J2 = 0.8, we have found the zero temperature param- 
eters m = 0.83, g x = 0.96, f y = 0.91, g x + y = 1.07, and 
Aj = 10.54J2. Note that, at T = 0, our calculation is 
consistent with that of RefAi. More details regarding the 
renormalizcd p and Aj obtained from a sigma model cal- 
culation will be discussed in Sec. V. Above T ct0 , the near- 
est neighbor bond correlations vanish, and two gaps be- 
come equal, Ai = A2 = ^/~p(8J2g x + y — p)- As the Ca v 
symmetry is restored above T a o, the velocity anisotropy 



disappears and v\ 3 



Vu 



V2i 



+y L 



IV. DYNAMIC STRUCTURE FACTOR 

The dynamic structure factor is calculated in the mod- 
ified spin wave theory through the average of the longi- 
tudinal and transverse spin structure factors. It is ex- 
pressed as 



1 



S ^u) = E [ c °sh(20 M 

k s,s— ±1 



2ft< 



ss\ 



x5(u - se k+q - s£k)«-k+q n k 



(21) 



?1 k = ?l k + 1 and = rik- 
Consider first uj -C T, and low temperatures T <C 
Aj. The dominant contribution to <S(q, uS) comes from 
the vicinity of the (n, 0) wave vector. In the limit 



ix\ « K 



T/v x and q y < A 



T/vy, we 
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FIG. 3: (Color online) The sharpening of the dynamic struc- 
ture factor around (n, 0) with decreasing temperature for 
u) = 0.3J 2 , S = 1, and J1/J2 = 0.8. 



can analyticall y 31 ' 33 ' 34 calculate <S(q, u>), which satisfies 
a dynamic scaling relation 

S(ir - q x ,q y ,w) = tS (tt - q x , q y )®(z, cjt), (22) 

where So(tt — q x ,q y ) is the equal time structure factor, 
and t = A^ 1 is the scaling time. So also satisfies a 
scaling form S {tt - q x ,q y ) = £ x £ y /(4ir\ y )A(z), where 

z = [{ixi. 11 — Qx) 2 + 6y9ij] 1 ^ 2 /2- The scaling functions are 
given by 



1 ( Q(x 2 - y 2 )- 

2k{x)\yWx 2 + {x 2 -y 2 ) 2 \ 'tt 



x arctan 



2 2 
x — y z 



x 2 + (x 2 — y 2 ) 2 



+ Q(y 2 -x 2 -l) 



A(z) 



logjg + Vj + z 2 

Zy/l + Z 2 



(23) 



When z -> 0, A(z) -> 1, and for z > 1, A(z) -> 
log(z)/z 2 . The second limit corresponds to momen- 
tum scales between inverse correlation length and inverse 
thermal length, where the system appears to have long 
range order (Goldstone mode behavior). The results for 
Sq are in agreement with one loop scaling results of a 
quantum nonlinear sigma mode l 35 ' 36 . 

A number of features follow from Eqs. (|22I23|) . As a 
function of energy for a fixed q with z 3> 1, <S(q, w) 
has a broad peak around ui ~ z/t. As a function of 
q for a fixed w, iS(q, ui) sharpens as temperature is re- 
duced reflecting the increase of correlation length; this 
is also seen from the results of direct numerical calcula- 
tions (Fig. In the numerical calculations of <S(q, ui) in 
Eq. (|2ip , a Lorentzian broadening of the delta functions 
has been employed, and consequently the gap between 
lot < z and lot > \/ z 2 + 1 is not observed in Fig. [3] 
but is instead left as shoulders. The processes beyond 
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FIG. 4: (Color online) Distribution of the dynamic structure 
factor in the momentum space for different temperatures and 
energies. The temperatures and frequencies corresponding 
to panels (a)-(d) are respectively given by (a) T/J2 = 0.5, 
L0/J2 = 2.0, (b) T/J 2 = 2.1, lu/J 2 = 2.0, (c) T/J 2 = 2.1, 
u)/J 2 = 3.0, (d) T/J 2 = 2.1, L0/J2 = 4.5. 



the modified spin wave theory are expected to smear the 
two-peak structure and also modify the scaling time r to 
the phase coherence time ~ (Aj/T) 1 / 2 /Ai^^. 



Beyond the u> <C T limit, we focus on the distribution 
of spectral weight in momentum space. Figs. 4(a) and 
4(b) illustrate the behavior at low energies. Provided 
T < T CT o, the anisotropy of the correlation lengths gives 
rise to an elliptic feature centered around (it, 0). The 
overall size of the ellipses is reduced as the temperature 
is decreased, reflecting increasing correlation lengths. On 
the other hand, the cllipticity has only weak temperature 
dependence; the ratio of two correlation lengths is almost 
unaffected by temperature variations for T <C T a o, due 
to the weak temperature dependence of the velocity ratio 

Vlx/Vly. 



With increasing energy, the evolution of the spectral 
weight distribution is illustrated in Figs. 4(b)-4(d). At 
intermediate energies, when u is comparable to the peak 
energy in the dispersion ek (see Fig. 2), there are fea- 
tures near ((1 ± 1/2)tt, 0), whose spectral weight is rel- 
atively small at the temperature shown in Fig. 4(c) but 
will increase with lowering temperature. The most visible 
spectral feature, however, is associated with the expand- 
ing ellipses surrounding (±7r, 0) and (0, ±7r), as is clearly 
seen in the high-energy spectrum shown in Fig. 4(d). 



V. THE ROLE OF ITINERANT ELECTRONS 
AND GINZBURG-LANDAU CONSIDERATIONS 

A. Anti-ferromagnetic fluctuations 

The description of the iron pnictides in terms of bad 
metals invokes quasi-localized moments coupled to itin- 
erant electrons whose spectral weight depends on the 
proximity of the system to the Mott transition^. For 
the parent compounds, the low-energy spin dynamics can 
be described in terms of a Ginzburg-Landau functional^ 
S = S-2 + S4 + . . . , where 



S 2 = J dqdu;[(r + wAq + cq 
+v{ql - ?y)m • m'], 



7 M)(m 2 + rrO 
(24) 



where m and m' are 0(3) vectors respectively for the 
magnetizations of the two decoupled sublattices, q x and 
q y are measured with respect to (±7r, 0) or (0, ±7r), w < 1 
is the coherent fraction of the single-electron spectral 
weight, and 7 is the strength of spin damping caused 
by the coupling to the itinerant electrons. £4 contains 
not only terms of the form m 4 , m' 4 and m 2 m' 2 , but 
also an order-from-disorder term (m • m') 2 with a nega- 
tive coefficient^. Eq. (|24| implies that elliptic features 
will occur in the dynamical responses even in the regime 
where the Ising order is not static but fluctuating and 
short-ranged; the primary role of the itinerant electrons, 
beyond shifting r through the positive wAq term, is to 
provide damping effects to such features. 

Well below the mean-field Ising transition tempera- 
ture, the thermal fluctuations of the Ising order param- 
eter a = ±(m ■ m )/|m||m | in the effective action of 
Eq. (|2~4"]) can be ignored. The choice of a = ± respec- 
tively correspond to short range (7r,0) or (0,7r) order. 
For short range (ir, 0) order, m — m becomes gapped 
and we find that order parameter dynamics can be ap- 
proximately determined in terms of a single 0(3) order 
parameter field M = m + m . The effective action for 
this field at quadratic order is given by 



So 



T 



J 1 



■j A 



Q 



2 2 
Qx v x 



2 2 
7 V 

h y 



7M] 



where v^ y = (c ± v/4), and u>i = 2irTl is the Matsubara 
frequency. With further assumption of small amplitude 
fluctuations, we can write M = Mon, where Mq is the 
constant amplitude, and n is the unit vector field. Thus 
low energy dynamics is now determined by a damped, 
anisotropic nonlinear sigma model. We consider the fol- 
lowing damped nonlinear sigma model action 



T 

2w? 



d 2 q 



V 



7M] |n(q,Wi) 



(26) 
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In writing the above equation we have rescaled 
yjv y /v x q x -> q x , and y/v x /v y q y -> g y , and w = y/Vx~Vy~, 
to write the action in spatially isotropic form, and g = 
v I p = v -1 xZ is the coupling constant with dimension of 
length. The scaling behavior of the correlation length in 
the quantum disorder phase and quantum critical regime 
for this damped nonlinear sigma model has been analyzed 
in Ref. l44l Here we will only consider the thermally dis- 
ordered or renormalized classical regime. In the large 
TV limit, the gap in the excitation spectrum A can be 
determined from the saddle point equation 



to obtain the gap A at large Landau damping, 



TV 



d 2 q 



1 



^ l A (2tt) 2 v*q 2 + uf + 7 M + A 2 



1 

vg 



(27) 



where A ~ n/a is the momentum cutoff. The Matsub- 
ara sum can be performed in terms of digamma func- 
tions, and after the momentum integration the left hand 
side can be expressed in terms of the logarithm of the 
gamma function. Here we consider two extreme limits of 
7/(2ttT) < 1 and j/(2nT) > 1. 

In the limit 'j/(2ttT) C 1, wc obtain z = 1 nonlinear 
sigma model result 



. . A , , »A / 2nv 

sinn — = sinn — cxp 

2T 2T P V 9 T 



(28) 



In the limit of small temperatures, such that vA T, 
and A <C T, wc obtain the result for small 7 limit, 



2nv (I 1 

A = Texp -— 

T \g g c \ 



A = v A exp 



2nv 
~T~ 



(1/g - l/g C 2) = wAexp 



¥) 

(31) 

with l/g c2 - (wA 2 )(21og2 - 1)/(4tt 2 7). Notice that cor- 
rect renormalized classical behavior of the correlation 
length for large Landau damping is found from the sad- 
dle point equation. For 27rp <Twc find z = 2 quantum 
critical behavior A ~ ^fhrfl \ augmented by logarithmic 
corrections. From the expressions of g c i, g c i we find that 
the stiffness for the z = 2 case is smaller than the z = 1 
case. This reflects the role of Landau damping. 

To summarize, in the limits of both small and large 
Landau damping, the correlation length has an exponen- 
tial temperature dependence in the renormalized classical 
regime. This will be the basis of our fitting the correla- 
tion length, which is described in Appendix A. 

If we consider the effects of the inter-layer antiferro- 
magnetic exchange coupling J z in addition to the J\ — J2 
model by using modified spin wave theory (see Appendix 
C), we obtain a finite mean field anti- ferromagnetic tran- 
sition temperature T/vo- Within the Ginzburg-Landau 
framework this corresponds to setting r(T) = 0. The 
fcrmion contribution wAq being positive, will decrease 
the transition temperature from the mean field value T/vo 
to a smaller value Tpf. However there will be significant 
amount of three dimensional antifcrromagnctic fluctua- 
tions up to the mean field Neel temperature T/vo- Above 
T/vo the magnetic fluctuations are essentially two dimen- 
sional. 



where g c \ = At:/ A is the coupling strength for zero tem- 
perature 2 = 1 quantum critical point, and p is the 
renormalized spin stiffness constant. From this expres- 
sion wc find £ = v/T cxp( 2 ^£) in the renormalized clas- 
sical regime described by T <C 2np. For 2np ^> T, 
one obtains z = 1 quantum critical behavior £ ~ v/T. 
If we go beyond the N — > 00 limit, or perform a two 
loop renormalization group calculation in the renormal- 
ized classical regime, we will find the correct classical 
result £ ~ exp(2f£)2i^ . 

For j/(2ttT) ^> 1, the physical properties are gov- 
erned by a z = 2 nonlinear sigma model. The fre- 
quency sum is performed after imposing a frequency cut- 
off ui c = ?j 2 A 2 /7, and after performing the momentum 
integration we obtain 



2nv , vA . 
= log — + log T 1 

gT A 



2i> 2 A 2 A 2 \ 
271-7T + 271-7T ) 



v 2 A 2 A 2 \ / 

i+ 2^ + 2^) +iogr v 1+ ^; 



(30) 



Now in the limit 2v 2 A 2 /(2tt7T) > 1 and A 2 /(27r 7 T) < 
1, we can use the asymptotic behavior of the log L(l + x) 



B. Ising fluctuations 

Since the Ising order parameter breaks Ci v symme- 
try, and in particular corresponds to B\ g representa- 
tion of the tetragonal lattice, it will couple to all the 
singlet fermion bilinears, which correspond to B\ g rep- 
resentation. Without the loss of generality if we con- 
sider a two orbital model of fcrmions including only 
dxz and d yz orbitals, the Ising order parameter a will 
couple to (cosfc.,. - cosfc 9 )* ks *k s , *j cs r 3 *k S) (cosk x + 
cos ky^^T^^s, and cos k x cos ky^^T^^ etc -' wnere 



( c t 

V xz, 



ks> yz,ks 



) describes the orbital and spin de- 
pendent fcrmion creation operators, and the Pauli matrix 
T3 acts on the orbital basis. Among the various B\ g bilin- 
ears, the conventional nematic order parameter and the 
ferro-orbital order parameter respectively correspond to 
(cosk x — cosfc !/ )^I'j cs ^I'ks and V&^ra^ks- Notice that wc 
can couple other d orbitals, following the same symme- 
try based criterion. When we integrate out the itinerant 
fcrmions, the contributions to the Ising order parameter 
a will arise from generalized _B lg particle-hole suscepti- 
bilities, and the quadratic part of the low energy action 
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for a will have the form 



wAn 



Ml 
q 



Wi,<»i)\ 2 

(32) 



In the above equation 7 CT is the Landau damping 
strength, and r a is the mass term arising from the lo- 
calized model, and wAq > is fcrmion contribution to 
the Ising mass. This fcrmionic contribution will suppress 
the Ising transition temperature from its mean held value 
T a0 to T a . But, the correlation length of the Ising order 
parameter will remain appreciable up to the mean held 
temperature T a0 . Since Ising transition occurs due to in 
plane magnetic fluctuations, consideration of inter-layer 
coupling does not signiheantly modify the Ising correla- 
tions. 

When we consider the magnetic and Ising order param- 
eter fluctuations on the same footing, further changes in 
the transition temperatures will arise from the self inter- 
action of a, m, m , and their mutual interaction erm-m . 
The interplay of Ising and magnetic order parameters, 
and their self-interactions arc crucial to determining if 
there will be a concomitant first order transition or two 
separate second order phase transitions. Despite the sup- 
pression of actual transition temperatures and the possi- 
ble complexity regarding the actual nature of the transi- 
tions, we still expect that the correlation lengths of the 
magnetic and the Ising order parameters will remain siz- 
able up to their respective mean-field transition temper- 
atures. 



VI. IMPLICATIONS FOR IRON PNICTIDES 

Our detailed theoretical studies provide the basis to 
understand the anisotropic spin responses that have 
been observed in the paramagnetic phase of the par- 
ent iron pnictides CaFc2As2^. These observations, 
made at temperatures above the first order antiferro- 
magnetic/structural transition, can be understood if the 
transition temperature is assumed to be considerably 
lower than the mean-field Ising transition temperature 
by the effects of fluctuations and coupling to phonons. 
To compare our theoretical results with the experiments 
of Rcf. |27| we have fitted the low frequency experimen- 
tal data with the dynamic structure factor calculated 
within the saddle point approximation of an anisotropic, 
damped nonlinear sigma model, which follows from the 
action of Eq. ([24]) . Within the saddle point approxima- 
tion the imaginary part of the staggered susceptibility is 
given by 



X (q-Q,w) = 



7 2 o> 2 + (cj 2 - v 2 (q x - it) 2 — v 2 q 2 — A 2 ) 2 

(33) 

The velocities of the effective model arc taken from the 
modihed spin wave calculations. The details of our pro- 
cedure are provided in Appendix A. 
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FIG. 5: Panels (a) and (b) respectively demonstrate <S(q — 
Q,oj = 12meV) at T = 1&0K obtained from our theory and 
data of Ref. [27]. Panels (c) and (d) respectively demonstrate 
5(q — Q,^ = 39meV) at T — 180-K" obtained from our theory 
and data of Ref. [27|. We have used J1/J2 = 0.552£, J 2 = 
9.8meV and 7 = 47meV. To facilitate the comparison with 
experimental result, we have plotted here in the Brillouin zone 
corresponding to the two-Fe unit cell instead of that for the 
one-Fe unit cell used in the rest of the paper. 



The comparison of our results with that of Rcf. [27| 
are shown in Fig. 5. The calculated elliptic features of 
<S(q — Q, ui) (Fig. 5(a)) is compatible with that seen ex- 
perimentally (Fig. 5(b)) at low frequencies. This con- 
tinues to be the case at higher frequencies, as shown in 
Fig. 5(c) and Fig. 5(d). The experimental results in the 
paramagnetic phase are consistent with our conclusions 
that as temperature is lowered, the peaks in the momen- 
tum space sharpen but the cllipticity is only weakly af- 
fected. Our estimated values of exchange constants are 
consistent with that of Ref. [I?]- When u> is smaller than 
the excitation gap A, the dynamic structure factor is 
peaked at q = Q. For lu > A, the intensity peak gets 
shifted to |q — Q| = Vw 2 — A 2 /v as shown in Fig. 5(c), 
and the ui 2 term in the dynamics is important to capture 
this feature also observed in the experiment as shown in 
Fig. 5(d). 

Inter-layer magnetic couplings in the parent iron ar- 
senides vary considerably amon g th e materials, but are 
always relatively weak. In Ref. [27|, the inter-layer cou- 
pling J z in paramagnetic CaFe2As2 was shown to be 
very weak, with J z /J2 = 0.1, being smaller than its 
counterpart in the magnetically ordered phase at low 
temperature. Consideration of such a weak inter-layer 
coupling does not appreciably change the estimated ex- 
change constants and the in-plane spin dynamics. An 
estimation of the spin stiffness constant using a renor- 
malized classical approximation for the correlation length 
shows that both fermion induced moment reduction, and 
Landau damping can sufficiently rcnormalize the stiff- 
ness constant (see Appendix A). In Appendix B we have 
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considered the effects of the weak inter-layer exchange 
coupling J z using the modified spin wave theory. For 
Jz/Ja = 0.1 the mean field Neel temperature Tjvo and the 
mean field Ising transition temperature T a @ become very 
close. However as we have discussed in Sec. V, despite 
the suppression of the actual transition temperature due 
to various fluctuation mechanisms, the magnetic and the 
Ising correlation lengths remain sizable up to the mean- 
field transition temperatures. In the temperature regime 
Ttv < T < Ttvo, there are three-dimensional antifcrro- 
magnetic fluctuations. However if we consider the ratio 
of the in-plane and inter-plane correlation lengths (mea- 
sured in units of corresponding lattice spacing), we find 
6s/£z ~ (Jz/{2J2 + Ji)) 1 / 2 . This ratio is of course mate- 
rial dependent. For weak inter-layer coupling of Ref. l27l . 
this ratio is ~ 0.2, and magnetic fluctuations are indeed 
quasi-two dimensional. 

Finally our discussion regarding the effect of itinerant 
electrons is most pertinent to the parent systems, but 
is consistent with the experimental observation of simi- 
lar low-energy anisotropic responses in the carrier-doped 
iron pnictides^Zr— . 




150 200 250 300 



T(K) 

FIG. 6: (Color online) Comparison between the correlation 
lengths £ = yfalw extracted from our fitting and that of 
Ref. [U The fitted spin stiffness constants are shown in the 
inset. 



VII. SUMMARY AND CONCLUSIONS 

We have addressed the spin dynamics in the paramag- 
netic phase of a two dimensional J\ — J 2 antiferromagnet 
on a square lattice at a finite temperature, using mod- 
ified spin wave theory. Within the modified spin wave 
theory we have identified a mean field Ising transition 
temperature T aQ , below which the C^ v symmetry of the 
square lattice is spontaneously broken. In the Ising or- 
dered phase the system demonstrates short range (tt, 0) 
or (0, tt) antiferromagnetic order. In order to systemat- 
ically understand the finite temperature spin-dynamics 
in the paramagnetic phase of iron pnictides, we have de- 
scribed the fermionic contributions and self-interaction 
effects of the order parameter fields within a Ginzburg- 
Landau framework. We have found that the fermion con- 
tribution and the self-interaction effects can considerably 
decrease the Neel and the Ising transition temperatures 
from their corresponding mean field values. However the 
correlation lengths of the magnetic and Ising order pa- 
rameters can remain appreciable up to the mean field 
transition temperatures. Based on this assumption, we 
have fitted the experimental data of Ref. H3, using our 
theoretical results. The calculated anisotropic features 
of the spin response are compatible with experiments for 
different frequencies. 

Finally, our calculations of the spin fluctuations at 
high energies should help understand future experiments. 
High-energy spin spectrum at the low-temperature or- 
dered state of CaFe 2 As2^ has already provided valuable 
information on the x-y anisotropy of the exchange in- 
teractions. Similar experiments have recently been re- 
ported in BaFe 2 As2^ and SrFe 2 As2^, including at tem- 
peratures just above the Neel transition where strong or- 



bital anisotropy has develope d 1 — . It will be instructive 
to experimentally map out the high-energy spectrum at 
higher temperatures in the paramagnetic phase. 
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Appendix A: Procedure of comparing with the 
experimental data 

We analyze the experimental data of Ref. by us- 
ing the imaginary part of the staggered susceptibility 
%(q — Q, w), which is calculated within the saddle point 
approximation for Eq. 1261 From the saddle point calcu- 
lation we find 

x"(q-Q w ) = glVj 

7 2 cj 2 + (oj 2 — v 2 (q x — tt) 2 — v 2 q 2 — A 2 ) 2 

(Al) 

At q = Q, we have 

X (0, w) = 2 A2 , 2 (A2) 

where x± is the uniform transverse susceptibility. We cal- 
culate the velocities v x and v y using Eq. [T3] and Eq. [TJ1 
and for the exchange constants we choose J 2 = lOmeV 
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and Ji = 0.55J2 = 5.5meV, as determined by Diallo 
et al£^. By fitting the experimental data we determine 
the temperature independent Landau damping strength 
7 and the temperature dependent gap A. By fitting the 
data for \ (0,w) at T = 1807^, with the formula from 
Eq. IA21 we find the Landau damping strength 7 and 
the gap A at 180A'. At low frequencies, Eq. (|A2[) can 
be further approximated by a Lorentzian with a width 
Tt ~ A 2 / ' \J^j 2 — 2A 2 . At the relatively low temperature 
of f 80K, the Lorentzian form is a good fit to Eq. (|A2|) up 
to frequencies of about 40 meV, and Ft = 7meV. At the 
high temperature of 300-fT, the Lorentzian form, which 
becomes a poorer fit to Eq. (|A2[) over the same frequency 
range, yields Ft = 4£meV. The definition of the energy 
line- width as Tt is the same notation as used in Ref. 1271 . 
but the constant 7 used in Ref. [2?] is not the conventional 
Landau damping strength and has a different meaning 
from ours. Our estimation is 7 = AlmeV . For the avail- 
able data at nonzero q — Q at different temperatures, we 
use the value of 7 so determined, and find the A at dif- 
ferent temperatures. Using the values of v x , v y and A, 
we find the correlation length £ x and £ y . The comparison 
of our theoretically calculated dynamic structure factor 
with the fitted parameter values, and the experimental 
results at low frequency 12meV are shown in Fig. 5(a) 
and Fig. 5(b) of the main text. Even at higher energy 
lj = 39meV our results for the dynamic structure factor 
are in reasonable agreement with experimental data, and 
the comparison for this frequency is shown in Fig. 5(c) 
and Fig. 5(d) of the main text. The consideration of the 
w 2 term in the effective action, leads to an interesting 
feature of the dynamic structure factor. For low frequen- 
cies such that oj < A, x (q — Q, w is peaked at q = Q. 
But at higher frequencies such that to > A, the intensity 
peak occurs away from the antiferromagnetic wave- vector 
and its location is determined by |q— Q| = Vw 2 — A 2 /v. 
This shift in the intensity peak can clearly seen by com- 
paring Fig. 5(a) and Fig. 5(c). Similar shift in the inten- 
sity peak can also be seen in the experimental results by 
comparing Fig. 5(b) and Fig. 5(d). For the correlation 
length we have compared our and experimental results 
in Fig. 6 by plotting the temperature dependence of the 
geometric mean of £ x and t; y . By fitting the correlation 
length with the renormalized classical formula, we have 
obtained an estimation for the stiffness constant. The 
fermion induced reduction of the magnetic moment Mo, 
and Landau damping are found to significantly reduce 
the stiffness constant in comparison to a pure J\ — J% 
model. 



we first incorporate the three-dimensional effects in our 
modified spin wave theory calculations. For simplicity 
we assume the sublattice angle <p = t- The modification 
to our discussion in Sec. I comes through an additional 
inter-layer antiferromagnetic bond correlation parameter 
g z . The Ising transition will be determined by the van- 
ishing of in plane nearest-neighbor bond correlations g x 
and f y . In the presence of J z , there is a finite, mean- 
field antiferromagnetic transition temperature Tjvo, cor- 
responding to Bose condensation of a's. The expression 
for total energy in Eq. [5] changes into 



E = - 



E (s+ 5-/(0) 



JiN 



2 \ 2 

Ol— ±X 



JiN 



E (s+l-m+fy 



02=±y 



E (s + \-f{Q) + g x+y 



2 ^ \ 2 

S3 — ±£±y 



J Z N 



2 ^ V 2 



(Bl) 



and the expressions for A k and B k are modified according 
to 



A k — 2Jig x C x ^ + 4J 2 g x + y C x+y ^ + 2J z g z C z ^ (B2) 
B k = 2Ji(g x - f y ) + 2Jif y CyM + 4J 2 g x+y - [i + 2J z g z , 

(B3) 

where C Zj k = cosfc z c, and c is the inter-layer separation. 
After accounting for the possibility of a finite staggered 
magnetization below TVo, the mean- field equations are 
given by 



f y = m + ~r ( " k + 5 ) C v-^ 



N V £ k 
k 

1 X - Ay 



(B4) 



A k ( 1 \ „ 
+ 2.^ — yik + ^ J c ",k, a = x, x + y, z 



5 +5 = m ° + ^E — 



, rikH — 
N *r< e k V k 2 

k 



(B5) 
(B6) 



Appendix B: Effects of inter-layer exchange coupling 

The quasi-2D nature of the spin dynamics was clearly 
shown in Ref. [2?]. To explain the observed (n, 0, tt) an- 
tiferromagnetic order, an inter-layer antiferromagnetic 
coupling J z was assumed and J z was estimated to be 
~ O.1J2. To assess the effects of J z on the spin dynamics 



For S = 1, J1/J2 = 0.8, c = a. the dependence of Tn Q 
and T a Q on J z / J2 are shown in Fig. 7(a). The temper- 
ature dependence of the mean-field bond parameters for 
•hi J2 = 0.1 are shown in Fig. 7(b). With increasing J z , 
the Neel temperature gradually increases and asymptot- 
ically approaches T a Q. Since the mean-field Ising tran- 
sition is a consequence of the two-dimensional magnetic 
fluctuations, T CT o is not modified by the finite inter-layer 
coupling J z . For J z / J2 = 0.1, and and J2 ~ lOmeV we 
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vj 2 




FIG. 7: Panel (a) shows the comparison between the mean- 
field Neel temperature Tjvo and mean-field Ising transition 
temperature T CT o, as a function of the inter-planar coupling J~, 
for JxjJi = 0.8, and 5 = 1. Panel (b) shows the temperature 
dependence of different mean-field parameters for J z /J2 = 
0.1. 



obtain Tjvo ~ T a Q ~ 2A0K, which is much higher than 
the actual Neel and structural transition temperature. 
Therefore the fluctuating anisotropy effects will be im- 
portant over a wide range of temperature, and the finite 



J z does not change this conclusion. 

Below T CT o, by expanding the dispersion around Q 
(-7T, 0, 7r), we obtain 



A 



[vKir - k x f + v 2 y k 2 y + v 2 z (n - k z f + A 2 ] 5 (B7) 

i 

[-H(8 hg x +y + ^Ji9x + 4 J z g z - fj,)] 2 , 

/i = 0,for T <T N (B8) 

(4J 2 g x+v + UiQxji^hgx+y + ^J\g x 



+2 J z g z ) 



{^hgx+y + 2Jig x + 2J 2 g z )(4J 2 3 a;+ y - 2Ji/ y ) 



B9) 



c p^^+y + 2^^ + 2J z g z )2J z g z }' : 



-2^010) 
(BH) 



We further notice that the velocities are well approxi- 
mated by 



2Sa(J 1 + 2J 2 )W1 



Ji + 2J 2 



2J 2 - Ji 
2J 2 + Ji 



a V 2J 2 + Ji 



(B12) 
(B13) 
(B14) 



and even in the presence of finite J z , the ratio v y /v x re- 
mains unchanged. For J z /J 2 = 0.1 and c/a rj 3.0262i we 
obtain v z /v x ~ 0.6, and this leads to smaller inter-planar 
correlation length (£ z < i/^Cy)- I R our comparison 
with experiments we have looked at the data that corre- 
sponds to in-plane dynamics, i.e., q— Q = (q x — tt, q y , 0), 
and consequently all the formulas remain unaffected. We 
also note that the effects of inter-planar coupling inside 
the magnetically ordered phase have been considered in 
Refs4^~— using similar technique. However our results 
are derived for the paramagnetic phase, which are essen- 
tially different from those described in Refs.— _ — . 
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